Josephson effect in ballistic graphene 
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We solve the Dirac-Bogoliubov-De-Gennes equation in an impurity-free superconductor-normal- 
superconductor (SNS) junction, to determine the maximal supercurrent I c that can flow through 
an undoped strip of graphene with heavily doped superconducting electrodes. The result I c ~ 
(W/L)eAo/h is determined by the superconducting gap Ao and by the aspect ratio of the junction 
(length L, small relative to the width W and to the superconducting coherence length). Moving 
away from the Dirac point of zero doping, we recover the usual ballistic result I c — (W/\F)eAo/h, 
in which the Fermi wave length \f takes over from L. The product 7 c J?n — Ao/e of critical current 
and normal-state resistance retains its universal value (up to a numerical prefactor) on approaching 
the Dirac point. 

PACS numbers: 74.45.+C, 74.50.+r, 73.23.Ad, 74.78.Na 



While the Josephson effect was originally discovered 
in a tunnel junction^ any weak link between two su- 
perconductors can support a dissipationless current in 
equilibrium^ The current I((f) varies periodically with 
the phase difference <f> of the pair potential in the two 
superconductors, reaching a maximum I c (the critical 
current) which is characteristic of the strength of the 
link. A measure of the coupling strength is the resis- 
tance i?N of the junction when the superconductors are 
in the normal state. The product I c i?N increases as 
the separation L of the two superconductors becomes 
smaller and smaller, until it saturates at a value of or- 
der Ao/e, determined only by the excitation gap Ao 
in the superconductors — but independent of the cou- 
pling strength. This phenomenology has been well es- 
tablished in a variety of superconductor-normal-metal- 
superconductor (SNS) junctions 3 and forms the basis of 
operation of the Josephson-Field-Effect-Transistor4£ 

A new class of weak links has now become available for 
research^ in which the superconductors are coupled by a 
monoatomic layer of carbon (= graphene). The low- lying 
excitations in this material are described by a relativistic 
wave equation, the Dirac equation. They are massless, 
having a velocity v that is independent of energy, and 
gapless, occupying conduction and valence bands that 
touch at discrete points (= Dirac points) in reciprocal 
spaced Graphene thus provides a unique opportunity to 
explore the physics of the "relativistic Josephson effect" 
(which had remained unexplored in earlier work& on rel- 
ativistic effects in high-temperature and heavy-fermion 
superconductors). We address this problem here in the 
framework of the Dirac-Bogoliubov-De-Gennes (DBdG) 
equation of Ref. 0. 

The basic question that we seek to answer is what 
happens to the critical current as we approach the 
Dirac point of zero carrier concentration. Earlier 
theoriesa2»ilii2i have found that undoped graphene has 
a quantum-limited conductivity of order e 2 /h, in the ab- 
sence of any impurities or lattice defects. We find that 
the critical current is given, up to numerical coefficients 
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FIG. 1: Schematic of a graphene layer, partially covered by 
two superconducting electrodes (S). A dissipationless super- 
current flows in equilibrium through the normal region (N), 
depending on the phase difference between the two supercon- 
ductors. Separate gate electrodes (not shown) make it pos- 
sible to vary independently the carrier concentration in the 
normal and superconducting regions of the graphene layer. 

of order unity, by 

I c ~ ^ max(W/L, W/Xp), (1) 
n 

in the short-junction regime L <C W, £ (with £ = hv/A 
the superconducting coherence length, W the width of 
the junction, and the Fermi wave length in the nor- 
mal region). At the Dirac point Xp — > oo, so the criti- 
cal current reaches its minimal value of (eAo/7i) x W/L. 
Since the normal-state resistance has its maximal value 
i?N — {h/e 2 ) x L/W at the Dirac point, the I c Rm prod- 
uct remains of order Ao/e as the carrier concentration is 
reduced to zero. 

The system considered is shown schematically in Fig. 
n A layer of graphene in the x — y plane is covered by 
superconducting electrodes in the regions x < —L/2 and 
x > L/2. The normal region |x| < L/2 has electron and 
hole excitations described by the DBdG equation^ 

Here Hq = —ifiv{a x d x +<7 y d y ) is the Dirac Hamiltonian, 
e > is the excitation energy, and fi is the chemical 
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potential or Fermi energy in the normal region (measured 
with respect to the Dirac point, so that [i = corresponds 
to undoped graphene). The electron wave functions fy e 
and the hole wave functions ~^>h have opposite spin and 
valley indices, which are not written explicitly. (A four- 
fold degeneracy factor will be added in the final results.) 
The Pauli matrices Oi in Hq operate on the isospin index, 
which labels the two sublattices of the honeycomb lattice 
of carbon atoms. 

Andreev reflection at a normal-superconductor (NS) 
interface couples ^ e and ^/ h . This coupling may be de- 
scribed globally by a scattering matrix, as was done in 
Ref.0 to determine the conductance of an NS junction. 
Here we follow a different approach, more suited to deter- 
mine the energy spectrum (and therefrom the Josephson 
current). In this approach electrons and holes are cou- 
pled locally by means of a boundary condition on the 
wave function in the normal region. 

We consider the energy range e < A below the excita- 
tion gap Ao in the superconductor, where the spectrum 
is discrete. At a point r on the NS interface (with unit 
vector n pointing from N to S, perpendicular to the in- 
terface), the boundary condition takes the form 



* h (r) = Af* e (r), 
1 

A 



(3) 



M = — (e - i^\A\ 2 -e 2 n ■ a) = e -i*-if>n*r. (4) 



Here A = Aoe 4 * is the complex pair potential in S, 
c = ipx^y) is the vector of Pauli matrices, and (3 — 
arccos(e/A ) G (0,7r/2). 

The relation © follows from the DBdG equationjSiiii 
under three assumptions characterizing an "ideal" NS in- 
terface: I) The Fermi wave length \' F in S is sufficiently 
small that A' F <C £, Af, where \f = hv/fi is the Fermi 
wave length in N and £ = hv/A is the superconduct- 
ing coherence length; II) The interface is smooth and 
impurity free on the scale of £; III) There is no lattice 
mismatch at the NS interface, so the honeycomb lat- 
tice of graphene is unperturbed at the boundary. The 
absence of lattice mismatch might be satisfied by de- 
positing the superconductor on top of a heavily doped 
region of graphene. As in the case of a semiconductor 
two-dimensional electron gasji^i^ we expect that such 
an extended superconducting contact can be effectively 
described by a pair potential A in the x — y plane (even 
though graphene by itself is not superconducting) . 

The particle current density out of the normal region, 
given by 



jparticle = V$* e fl ■ CT* e - V^* h fl ■ <T#fc, 



(•5) 



should vanish for e < Ao, because subgap excitations 
decay over a length £ in S. (The possibility of a sub- 
gap excitation entering the superconductor at one point 
along the boundary and exiting at another point within 
a distance £ is excluded by assumption II.) By substi- 
tuting the boundary condition @ one indeed finds that 
Jparticle — 0, since M is a unitary matrix which commutes 
with n ■ er. 



In the SNS junction the normal region has two inter- 
faces with the superconductor, one at x = —L/2 (with 
superconducting phase $ = 0/2 and outward normal 
n = —x) and another at x = L/2 (with <!> = —0/2 
and n = x). The boundary condition (J3J at the points 
r± = (±L/2,y) thus takes the form 

* fc (r_) = tf(e)* e (r_), * h (r+) = C7- 1 (e)* e (r+X6) 



U(e) 



_ -i0/2-M/3cr x 



(3 = arccos(e/A ). 



(7) 



Since the wave vector k y parallel to the NS interface 
is conserved upon Andreev reflection, we may solve the 
problem for a given k y = q. The transfer matrix Ai(e, q) 
relates the states at the two ends of the normal region: 

* e (r+) = M(e, q)# e (r-), V h (r+) = M(-s, g)# fc (r_). 

(8) 

(For ease of notation, the ^-dependence will not be writ- 
ten explicitly in what follows.) The condition for a bound 
state (= Andreev level) in the SNS junction is that the 
transfer matrix for the round-trip from r to r + and back 
to r_ has an eigenvalue equal to unity. This condition 
can be written in the form of a determinant, 



Det [1 - M- 1 (e)U(e)M(-e)U(e)] = 0, 



(9) 



which we have to solve for £ as a function of q and <p. 

The electron transfer matrix M. (e) is readily obtained 
from the Dirac equation, 



M = ke lkLa *k, 



A = A" 1 = (2cosa)- 1 / 2 ' ' " ' 



ia/2 gia/2 
gia/2 _g-ra/2 



a e 



hvq 



k(e) — (hv)^ 1 (e + n) cosa(e). 



(10) 
(11) 

(12) 
(13) 



The angle a is the angle of incidence of the electron, and 
k is its longitudinal wave vector. 

Evaluation of the determinant © leads after some al- 
gebra to the quantization condition 
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cos 2/3 

sin 2(3 
(14) 



where we abbreviated a± = a(ie), 9± = k(±e)L. 

We introduce a finite width W to quantize the trans- 
verse wave vectors, q —* q n , n = 0, 1, 2, . . ., and denote by 
p n (e 7 4>) the density of states in mode n. The Josephson 
current at zero temperature is then given by 



oo 



(15) 



where the factor of 4 accounts for the two-fold spin and 
valley degeneracies. To be definite we take "infinite 
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mass" boundary conditions at y = 0, W, for which 11 
q n = (n+l/2)ir/W. (For W 3> L the choice of boundary 
conditions becomes irrelevant.) At the Fermi level, the 
lowest N(fi) — fiW/irhv modes are propagating (real k), 
while the higher modes are evanescent (imaginary k). 

We analyze the Josephson effect in the experimentally 
most relevant short-junction regime that the length L of 
the normal region is small relative to the superconduct- 
ing coherence length £. In terms of energy scales, this 
condition requires Ao <C hv/L. To leading order in the 
small parameter AqL/Hv we may substitute a± — > cv(0), 
6± — > k(Q)L in the quantization condition Q14|). The so- 
lution is a single bound state per mode, 



e n (4>) = A o ^l-r„sin 2 (0/2), (16) 

k 2 cos 2 (fc n L) + sin 2 (fc ra L) ' ^ 

with k n = [(fi/hv) 2 — g^] 1 / 2 - This expression for the An- 
dreev levels in terms of a normal-state transmission prob- 
ability r n has the usual form for a short SNS junction^ 
Comparison with Ref. [ll| shows that r„ is indeed the 
transmission probability for a ballistic strip of graphene 
between two heavily doped electrodes in the normal state 
(Ao = 0, X' F <C Af). The normal-state resistance i?N is 
thus given by 



-i 4e 2 ^ 



71=0 



(18) 



Substitution of p n (e,(f>) = S[e — e n {4>)} into Eq. ltT3|) 
jives the supercurrent due to the discrete spectrum, 



eA x - 
h 2^ 



r„ sin ( 



-[l-r n sin 2 (^/2)]V2 



(19) 



Contributions to the supercurrent from the continuous 
spectrum are smaller by a factor L/£ and may be ne- 
glected in the short-junction regime>i£ For L <C W the 
summation over n may be replaced by an integration. 
The resulting critical current I c and the I c i?N product 
are plotted as a function of fj. in Fig. 

The limiting behavior at the Dirac point (jx <C hv/L) 
for a short and wide normal region (L <C W, £) is 

I{4>) = ^ cosU/2) artanh [sm(<£/2)l, (20) 
n ttL 

eAn W 

/ c = 1.33 — ^— , 7 c i? N = 2.08A /e. (21) 

These results for ballistic graphene at the Dirac point are 
formally identical to those of a disordered normal metal 
(Fermi wave vector kp, mean free path Z) ) 16 ' 18 upon sub- 
stitution kpl — > 1. This correspondence is consistent 
with the finding of Ref. ITU that ballistic Dirac fermions 
have the same shot noise as diffusive nonrelativistic elec- 
trons. 
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FIG. 2: Critical current 7 C and I c Rn product of a ballistic 
Josephson junction (length L short compared to the width W 
and superconducting coherence length f), as a function of the 
Fermi energy \i in the normal region. The asymptotes 12 U 
and 1221 are indicated by dashed lines. 



In the opposite regime /i ^> hv/L we have instead (still 
for L < IF, £) the result 



1.22——^^—, 7 c -Rn 

ft 7mU 



2.44 A /e. (22) 



(We do not have a simple analytic expression for the 0- 
dependence in this regime.) The critical current l|22|) is 
about half the ideal ballistic valu o 16 ' 19 I c = 2NeAo/h, 
with N = fiW/irhv the number of propagating modes 
(per spin and valley). This reduction is due to the mis- 
match in Fermi wave length at the NS interfaces. Eqs. 
(|21|l and ll'L'l) together contain the scaling behavior (|TJ 
announced in the introduction. 

In conclusion, we have shown that a Josephson junc- 
tion in graphene can carry a nonzero supercurrent even 
if the Fermi level is tuned to the point of zero carrier 
concentration. At this Dirac point, the current-phase re- 
lationship has the same form as in a disordered normal 
metal — but without any impurity scattering. Instead 
of being independent of the length L of the junction, 
as expected for a short ballistic Josephson junction, the 
critical current I c at the Dirac point has the diffusion- 
like scaling oc 1/L. Since the normal-state resistance 
i?N oc L, the I c Rn product remains fixed at the su- 
perconducting gap (up to a numerical prefactor) as the 
Fermi level passes through the Dirac point. This unusual 
"quasi-diffusive" scaling of the Josephson effect in un- 
doped graphene should be observable in submicron scale 
junctions. 
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